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Renormalization-group (RG) flow equations have been derived for the generahzed sine-Gordon 
model (GSGM) and the Coulomb gas (CG) in d > 3 of dimensions by means of Wegner's and 
Houghton's, and by way of the real-space RG approaches. The UV scaling laws determined by 
the leading-order terms of the flow equations are in qualitative agreement for all dimensions d > 
3, independent of the dimensionality, and in sharp contrast to the special case d = 2. For the 
4-dimensional GSGM it is demonstrated explicitly (by numerical calculations), that the blocked 
potential tends to a constant effective potential in the infrared (IR) limit, satisfying the requirements 
of periodicity and convexity. The comparison of the RG flows for the three-dimensional GSGM, the 
CG, and the vortex-loop gas reveals a signiflcant dependence on the renormalization schemes and 
the approximations used. 

PACS numbers: ll.lO.Hi, ll.lO.Kk 

I. INTRODUCTION 

The renormalization of a two-dimensional one-component scalar field theory with a periodic self-interaction, and 
the two-dimensional generalized sine-Gordon model (GSGM) have been investigated for d = 2 dimensions in great 
detail over the last three decades [1111110,111113,111111111111111 

Here, we only mention some of the 

main results: 

• The blocked potential for the two-dimensional GSGM tends to a constant effective potential in the IR limit, 
independent of the field [ll,[ll|ll. 

• Following the procedure proposed in |16| , one can identically rewrite the partition function of the d-dimensional 
GSGM in the form of a d-dimensional static Coulomb gas, which we call here the equivalent Coulomb gas 
(EGG). For d = 2 it is well-known that the GSGM and the EGG, as well as the XY model describing classical 
planar spins belong to the same universality class _10j. The XY model has a dual theory in the sense of 17] that 
can be reformulated as a gas of topological excitations: the two-dimensional vortex gas (VG) for d = 2, and the 
three-dimensional vortex-loop gas (VLG) for d = 3 (see fl^ and Refs. therein). For dimensions d — 2, it was 
shown that the EGG and the VG can be transformed into one another by an appropriate duality transformation 
111 ^1 that inverts the temperature. Two-dimensional generalized models are well known where both 
the EGG and the VG are included as particular limiting cases [1 1 1 0| and are self-dual under the duality 
transformation. 

While the relations of the various models for two dimensions are well-established, those of the corresponding 
models for d = 3 are not completely settled (see e.g. 0). In particular, it is not proven that the VLG and the three- 
dimensional XY-model belong to the same class of universality. Therefore, the investigation of the renormalization of 
the GSGM and the EGG for d > 3 dimensions is essential for a further clarification of this issue. The purpose of the 
present work is threefold: 

1. to investigate the renormalization of the GSGM for dimensions d > 2 by the Wegner- Houghton renormalization- 
group (WH-RG) method 0| and to demonstrate numerically that the blocked potential for the 4-dimensional 
GSGM tends to a constant effective potential in the IR limit; and 

2. to perform an RG analysis for the EGG for dimensions d > 2 by the real-space renormalization group (RS-RG) 
method using the dilute gas approximation (DGA) and finally 

3. to compare the RG flows obtained for the GSGM, EGG and VLG [T9ll20l | in the three-dimensional case. 

Renormalization of the GSGM and that of the related models in d > 3 dimensions are of particular interest with 
regard to the numerous physical applications. In high-energy physics the haaron model of confinement is in need 
of a clarification of the IR behaviour of the GSGM in d = 4 dimensions. In low-temperature physics the superfluid 
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transition of He''^ belongs to the universality class of the three-dimensional XY model and it provides an exceptional 
opportunity for an experimental test of the RG predictions |2]| . In superconductivity the three-dimensional XY 
model (isotropic and anisotropic) was used to study the phase transitions in the three-dimensional flux-line lattice 
|22l l23j| . The three-dimensional XY model can also help us in understanding the flux motions in superconductors 
with extremely high anisotropy; these have been observed in recent experiments 24\. 

Due to the fact that the GSGM and the EGG are equivalent at the level of the partition function, their 
renormalization-group flows obtained by applying different RG methods are alternative approaches to tackle essentially 
the same problem. On the one hand, the WH-RG applied to the GSGM keeps the periodicity of the self-interaction 
potential in each inflnitesimal blocking step. There is a price to pay: one has to handle inflnitely many vertices in the 
process of searching the blocked action in the subspace of functionals periodic in the field variable. Furthermore, the 
requirements of periodicity and convexity |25| are laid upon the effective potential of the model which constrain it to 
a field-independent constant The WH-RG has the advantage that (i) it can be extrapolated to the IR limit and 
(ii) is able to treat the spinodal instability 26] . Its main disadvantage is that the truncation of the gradient expansion 
is unavoidable and that the implementation of the sharp momentum cut-off can induce unphysical features of the 
RG flow going beyond the local-potential approximation (LPA) [23|. By restricting ourselves to LPA, we ignore the 
wave-function renormalization, which is equivalent to neglecting the scale-dependence of the temperature beyond the 
tree-level and its physical consequences. 

On the other hand, the RS-RG applied to the EGG keeps intact the CG structure in each blocking step and 
therefore also the periodicity, albeit indirectly. It implements a sharp cut-off in real space and has been proven 
very adequate for applications to the two-dimensional Coulomb gas and to the three-dimensional vortex-loop gas 
[i3,|23- Nevertheless, it has the drawbacks that (i) no systematic approach is known to go beyond the dilute-gas 
approximation (DGA), and (ii) the sharp cut-off in real space, which is effectively a smooth momentum cut-off, is 
inadequate to treat the spinodal instability. Because our WH-RG approach will show that the latter always occurs, the 
interpretation of the extrapolation of our RS-RG results for the EGG to the IR regime is problematic. Nevertheless, 
it is interesting to compare the RS-RG results for the EGG in the UV regime to those for the gas of topological 
excitations (e.g. the VLG for d — 3). 

The phase transition of systems in the universality class of the two-dimensional Coulomb gas, known as the 
Kosterlitz-Thouless-Berezinskii phase transition, is a topological phase transition because there is no spatial long- 
range order for d = 2 dimensions at flnite temperature due to the Mermin- Wagner theorem |28j . For d> 3 dimensions 
the Mermin- Wagner theorem does not hold any more. Therefore, one expects a modification of the phase structure 
of the Coulomb gas in the number of dimensions d > 2 as compared to the case with d = 2. As shown below, the UV 
scaling laws in d > 2 for both the GSGM and the EGG share the property of a lacking Coleman fixed point, which 
by contrast is present for d = 2 0] . 

In Sect. ^]we introduce our notations and summarize RG flow equations obtained in the literature for the various 
models discussed in the next sections. The RG flow and the blocked potential for the 4-dimensional GSGM are 
determined numerically in Sect. IIIII by the WH-RG method. In Sect. IIVI the RG flow obtained by the RS-RG 
approach for the d-dimensional Coulomb gas is discussed. In Sect. 0the RG flow of the three-dimensional Coulomb 
gas (EGG) is compared to the RG flow obtained in for the three-dimensional GSGM, and to those for the three- 
dimensional VLG 20]. Finally, the results are summarized in Sect. IVII RS-RG flow equations for the d-dimensional 
Coulomb gas are derived in App. ^ 



II. GSGM AND RELATED MODELS 



A. GSGM and EGG 



First, we remind the reader that the partition function for the GSGM can be identically rewritten in the form of 
the partition function for a CG, that we shall call the equivalent CG (EGG). In its lattice regulated form the partition 
function of the GSGM is given as 



-^GSGM = / T^[(l>] cxp - 2^^'' X! ^^^^,y^y ~ "'^ cos(n/3^2;) 



(1) 



x^y X n—1 

where a denotes the lattice spacing, and d^^y = 2a~'^~'^{dSx.y — X])!=i Sv,x+n^,a) is the inverse propagator. Performing 
the transformation of the field variable tpx = P4'x, one finds 

,2d 



Z = / V[(p] exp 



^ ^ (fxdx^y^Py - X! X! COS(n(/3:r) 



X n—1 



(2) 
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This partition function can be identically rewritten in the form of the macrocanonical partition function for a CG 
applying the procedure proposed in [l6| : 

1. One expands the exponential factor of the integrand with the periodic potential in Taylor series, expresses 
cos{n(p) in terms of the exponential function, and introduces the integer valued variable, the charge a = 
±1, ±2, . . .; the result of these operations is 



\ X / iy=0 



E 



U(Ti • • • Ucr^ exp 



xi ,2:2 ,■ ■ ■ ,2;^ 



-,2d 



x,y X 



(3) 



where the charge density, px — '^^x,xj, has been introduced and u^cr = holds due to the Z{2) 

symmetry of the potential. 



2. Performing the Gaussian path integral, one obtains 

1^=0 



Uctj • • • Ua-^ ^xl',...,x!^ : 



(4a) 



XI ^X2^...^Xu 
(Ti,...,cr„5^0 



rfTi,...,<T„ 

XI,. ..,Xv 



exp 



1 1 " 

-TrlnG-1 - ^ a.a,G.„,^. 



(4b) 



with the (dimensionful) propagator Gx,y — a ^'^(d ^)x,y 
Including the self-interaction terms into the (dimensionless) fugacities, ' 



^a'^Ua <2xp[—^a'^(3^Gxi,xi] and replacing 



G^^a: by the value G{a) of the propagator at the cut-off distance a, the partition function takes the form 



where 



Jxi-x: = exp 



1 1 " 

--TrlnG-i - :^/3' E '^'^j'^---. 



Wo 



a'^Ug- exp 



(6a) 



(6b) 



with the dimensionless coupling 0"^ — a and the dimensionless propagator 



Gx.y ^ 



-C2 In 



a 

-d-2 



Cd 



I (d-2)|2;-j/| 



for d = 2 , 
77=2 for d > 3 . 



(7) 



Here, Cd = where ^Id denotes the entire solid angle in d dimensions, and a — 0{a). Its explicit value has been 
determined for d = 2 (c.f. 10]). For the sake of simplicity we shall take here a — a. The integration regions on the 
r.h.s. of Equation (|SJ) are discussed in App. ^ Eq. ||SJ) represents the grand-canonical partition function for a gas 
of static Coulomb charges. The temperature of this EGG is given as oc (/3^)~^, and the Fourier amplitudes u^r of the 
periodic self-interaction potential in the GSGM play the role of the fugacities of the various charges a — ±1, ±2, . . .; 
this mapping of the partition function of the GSGM to that of the EGG is exact for arbitrary dimensionality. 

In this work we shall check our RS-RG results for the EGG in higher dimensions against the RG-flow equations 
obtained in 7] with the same approach for d = 2. It is a little cumbersome to identify the various notations (see 
Tab. m for an explanation). Indeed, the two-dimensional EGG is an electric Goulomb gas in the sense of 0. Applying 
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the RG flow equations f?', Eqs. (2.22) and (2.23)] to the two-dimensional EGG with Z{2) symmetry (wcy = w-o-), one 
obtains: 

cr=l " V / a'^a,cr ^ ^ <t' = 1 

From the hnearized flow equations one recovers the Coleman fixed point 0] of the sine-Gordon model at (3'^ — Stt. 
Below, we also compare our RS-RG flow equations derived for the EGG for arbitrary dimensionality to Eq. ((HJ. 



m 


ECG 


\1\ 


VG 


1/3 


^7(2^) 


9 


2tiJ 


e 


(7 


m 


Q 


F(e,0) 




y(o,m) 


A, 



TABLE I: Identification of the various notations for the parameters of 
the two-dimensional ECG (first and second columns), as well as for the 
VG (third and fourth columns) according to Eqs. (2.1) and (2.2) in 
and our Eqs. (EJ, 0, and itTUl . 



B. Gas of topological excitations 

In general, dual transformations IT^ can be found and applied to almost any Abelian theory in any number of 
dimensions. The dual form of the original theory represents an intermediate step towards its third form when the 
partition function is expressed in terms of the topological excitations of the original degrees of freedom (see y^]), the 
gas of topological excitations. 

The d-dimensional XY-model is given by the partition function 



ZxY = / VIS^] - 1) exp 



L <V,z> 



(9) 



where the classical spin Sx is a unit- vector in the two-dimensional internal space; X]<y2> stands for the sum over 
pairs of nearest neighbour lattice sites. For the numbers d — 2 and d = 3 of dimensions the corresponding gases of 
topological excitations are well-known: 



1. The vortex-gas (VG) for d ^ 2 dimensions. The partition function of the XY-model can be rewritten as (see jlOj ) 



■ Aq^ CXp 



1 " ~ 



(10) 



This describes a system of static vortex charges interacting via the two-dimensional Goulomb potential 

Gxi,xj = ^n{\xi - Xj\/a) , 

where the fugacities of the vortices are defined as Ag — exp[—fiq^J/2]; q = ±1,±,2, . . . , and fi can be interpreted 
as the chemical potential of the vortex-antivortex pairs. The VG and the ECG are both Goulomb gases, but 
those of different physical origin. The vortices are the topological excitations of the free periodic scalar field. 
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whereas the Coulomb charges of the ECG are connected with the periodic self-interaction of the scalar field. 
Using the identification of various notations given in Table (last two columns) , one may interpret the VG as 
a magnetic Coulomb gas in the sense of Again, applying the RG flow equations 0, Eqs. (2.22) and (2.23)] 
for the VG with Z{2) symmetry {Aq — A^q), one finds 

°^=(^""^'^')^«+" ^ Aq,Aq^q,-jf+Vi)Aqf2Al,. (11) 

According to Table U the RG equations ((TT)l can be rewritten in terms of = 1/ J as 

«^=4.3f;,M^, J-^^(2-Z-AAq + n Y: Aq,Aq^q,-J^-^ + V^)Aqf:Al. (12) 

g=l \ P / q'^0,q ^ ^ q' = l 

2. The vortex-loop gas (VLG) for c? = 3 dimensions. For the three-dimensional XY model the dual theory has 
been constructed in the literature [l^l20l | and it turned out to be the gas of interacting vortex loops. Also the 
corresponding fiow equations have been derived for the parameters K (x J and y (the fugacity of the vortex 
loops) by RS-RG in the DGA: 

where L = L{a) corresponds to the long-range contribution. One may assume L to approach a constant in the 
IR limit, or to be weakly (logarithmically) divergent, e.g. L{a) = 1 in Ref. 19], L{a) = ln(a/ac) + 1 in Ref. [20l ]. 



III. RG FLOW OF GSGM 



A. Differential RG in momentum space 



Below we shall discuss the renormalization of the GSGM (Euclidean one-component scalar field theory with pe- 
riodic self-interaction) by means of the differential renormalization group approach in momentum space. In this 
approach, blocking transformations ^291 ] are realized by successive elimination of the field fluctuations according to 
their decreasing momentum in infinitesimal steps. The physical effects of the eliminated modes are encoded in the 
dependence of the coupling constants on the scale k above which the high-frequency modes are effectively eliminated. 
The elimination of the modes above the moving scale k is either complete as in Wegner's and Houghton's method 
(WH-RG) (sharp IR cut-off) or partial as in Polchinski's method js^] (smooth IR cut-off). Both of these RG 
methods keep intact the periodicity, which is the essential symmetry of theory. Nevertheless, choosing the WH-RG 
method is a necessity in our case: as shown below, a spinodal instability occurs at some critical momentum scale 
kc- Therefore, the elimination of the modes above and below this critical scale involves completely different physics. 
Renormalization of the couplings occurs due to small quantum fluctuations above the critical scale, whereas tree-level 
fluctuations of large amplitude dominate below that scale [2^. Thus, the application of the RG methods with smooth 
cut-off like that of Polchinski's method [s^, |^ 1131 is now unjustifled, because such methods treat the modes at all 
scales similarly. 

Above the critical scale kc the WH-RG method |3 provides an integro-differential equation for the blocked action 
Sk[(t)] as the functional of the scalar fleld 0. In order to solve this equation, one has to project it to a particular 
functional subspace. Therefore, one generally assumes that the blocked action contains only local interactions, expands 
it in powers of the gradient of the fleld, and truncates this expansion at a given order, for technical reasons. The 
drawback of the WH-RG method combined with the gradient expansion is that the truncation of the latter may 
introduce an unphysical high-frequency behaviour of the blocked action, since a sharp cut-off is used p^ . Here, we 
restrict ourselves to the leading order of the gradient expansion, the local-potential approximation (LPA). In order to 
get more insight in the RG evolution of the blocked action, it may turn out to be useful in the future to replace the 
gradient expansion with some more appropriate approximation scheme, like the cluster expansion proposed in [33| . 
In LPA the ansatz for the blocked action can be written as 



5fc = d 



(14) 
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with the blocked potential 14 (</>) for which the WH equation reads as follows |3J, |35|, |3a, ISTj, |33 : 

/ e + dlVki(t>)\ 

fc9fcl4(0)=-fc%ln( p-^j' (15) 

with ad = ^fld{2n)^'^ , where fl^ denotes the solid angle in d dimensions. In order to look for fixed-point solutions of 
(|15|l or to follow the scaling of an arbitrary potential, it is convenient to remove the trivial scaling of the dimcnsionful 
coupling constants and rewrite Eq. (|15|l as 

(^d-^^^d^ + kdk^Vk{4')= -^ittd ln(l + 9|Vfe(<^)) , (16) 

where the dimensionless quantities (f> = k ~ (p and Vk{4>) — k^'^Vk{4') have been introduced. 

Notice that the argument of the logarithm in Eqs. (|15l) and H16|l must be positive. If the argument vanishes or if it 
changes sign at a critical value fcc, given by k^ = —d'^Vk^{(j)), the WH equation ((T5|l loses its validity for k < kc- This 
is a consequence of the spinodal instability. At the scale kc zeros occur among the eigenvalues of the second derivative 
matrix of the blocked action, the restoring force vanishes for some modes of the fluctuations and the amplitudes of 
such fluctuations can grow to a large finite value and cause a tree-level renormalization "2^. Then, the tree-level 
blocking relation [2^ 



Vk-ski(f>) = rain 



1 



k'^p'^ + ^ I du Vk{(j> + 2pcos{TTu)) 
2 J_i 



p 

should be used instead of Eq. H15|l for the determination of the blocked action Vk-5k{<j>) at the lower scale k — 5k. 



(17) 



B. Global aspects of the renormalization of the GSGM 

As it is argued in Refs. fT5,'T4'| for d = 2 dimensions, the WH-RG procedure retains the periodicity of the potential. 
The argumentation holds for any number of dimensions: both Eqs. H16|) and H17|) which describe the infinitesimal 
blocking step 5k can be considered in the form Vk-Sk{4') — -^[14(0)] with some functional T. If the r.h.s. of this 
equation is periodic at the scale k due to Vk{4>) — Vfe(0-f A), then its l.h.s., i.e. the blocked potential Vk-Sk{4>)i should 
also be periodic with the same period A. This means that (i) the WH-RG in LPA keeps periodicity of the blocked 
potential, and (ii) the period A in the internal space does not depend on the scale. Note that these statements are 
only valid for the dimensionful potential as the local functional of the dimensionful field. 

Since the blocked potential Vk{4>) tends to the effective potential Vos{4>) in the IR limit fc 0, the effective 
potential should be periodic, as well. On the other hand, the effective potential is convex over the homogeneous 
field configurations. The double requirement of periodicity and convexity can only be fulfilled in a trivial manner, 
namely if Ves{4>) is a constant function. In order to demonstrate this general statement, we determine the blocked 
potential for the 4-dimensional GSGM numerically. 



C. RG equations 

Here, we apply the flow equations and p7(l for the GSGM. For the sake of simplicity, we consider periodic 
potentials with Z(2) symmetry, Vk{<j)) — Vfe(— 0). The Fourier-expansion of such potentials is given by 

OC CXD 

Vfe(0) = ^u„(fc)cos(n/3</)), Vfe(0) = ^{t„(fc)cos(n/30) (18) 

n— 1 n— 1 

where fi — 27r/A is independent of the scale fc, and the tilde indicates the corresponding dimensionless quantities. 

Inserting the Fourier-expansion l|18() into Eq. H16(l differentiated with respect to (/>, one obtains an equation with 
terms cither non-periodic or periodic in the field variable (p. The non-periodic terms occur due to the scale-dependence 
of (3{k) and the term —^-^(f)d^Vk{4>) in Eq. (|16|) . The obtained flow equation is satisfled if and only if the non-periodic 
and the periodic parts arc identical on both sides separately. The non-periodic parts yield: 

J2<ik)n (kdkm - ^^W) ^ sin(n/3 0) = 0. (19) 



7 



Using this relation and the equation for t he p eriodic terms, one gets the evolution equations for arbitrary number d 
of dimensions as given previously in Ref. jl2|: 



1. 



[d + kdk)vn{k) = aP''{k)n'vn{k) - -^\k) ^ A„,p(fc)(d + kdk)vp{k), (20) 



kdk/3^ik) = {d-2)(3^ik), (21) 

where Vn{k) = nun{k), and An^p{k) = \n — p\vn-p — {n + p)vn+p [note that in our previous definitions the analog 
of the first relation is Vn{k) = n/3u„(fc)]. Eq. follows from Eq. (fT^ and has the solution 

p^k)^ k^^-^P^. (22) 

This expresses that the periodicity of the blocked potential is kept, whereas the dimensionful period A — ^ remains 
independent of the scale. 

In order to solve the flow equations one has to express the derivatives kdkUn{k) occurring in Eq. (pn|l explicitly 
and truncate the Fourier-expansion of the potential at some mode n = N. This can be achieved by the following 
steps: 

1. Rewrite Eq. (|2()|l in terms of Un = ^w„: 

+ 2/3' ^",P ) kdkUpik) = -du„{k) + adp^{k)n^un{k) - -P^{k) ^ An,p{k)up{k) , (23) 
p=i ^ ^ p=i 

with An^p{k) = ^[{n - p)'^Un-p - (n +p)^u„+p]. 

2. Invert the matrix on the l.h.s. of Eq. (I23|) numerically in the subspace of the modes with n,p < N. 

For later comparison with the flow equations for the d-dimensional ECG, it is useful to perform the inversion of the 
matrix 1 + keeping the terms up to the order 0{A^), i.e. those up to 0{u^) : 

OO OO 

p—l p,r— 1 

In the region of the spinodal instability for k < kc, the RG flow should be determined with the help of the tree- level 
blocking relation (|17|l which can be rewritten for the dimensionful periodic potential as 

Vk-Sk{4>) = rnin 
p 

where Jg stands for a Bessel function, /3 =constant and is a homogeneous field configuration. It is worthwhile to 
note that Eq. H25|) is just the same for any number d of dimensions. Therefore, it can always be solved in the same 
manner as discussed in |l2j | for the case with d — 2. The position p — pk{<f) of the minimum should be determined 
separately for each constant field configuration at the given scale k. Then, the r.h.s. of Eq. (|25|l at the minimum 
as a function of (j) should be reexpanded in Fourier-modes, in order to determine the coupling constants Un{k — 5k) of 
the blocked potential at the lower scale k — 5k. The pure tree-level scale-dependence (|22ll should hold in both regions 
k > kc and k < k^ alike. 



OO 

E 



Un{k) cos{nP(j))Jo{2nPp) 



(25) 



D. UV scaling 

The trivial Gaussian fixed point = is the only one of the Wegner-Houghton RG equation 116|l in the subspace 
of periodic potentials. This is because the dimensionless coupling /3(fc) can only be independent of the scale for a 
periodic potential that vanishes identically, as Eq. H22() should hold for a constant /3. 

In order to determine the scaling operators and the critical exponents of the periodic field theory around the 
Gaussian fixed point for d> 3 dimensions, one can use the quasi-linear version of the flow equation H24|) : 



{d + kdk)un{k) = adi? ri^ Un{k). 



(26) 
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Together with Eq. (|22|l this has the analytic solution 



M„(fc) = Un{A) [ — ] exp 



r ad/32 (A)n2 



AJ d-2 



(27) 



where /3(A) and w„(A) are the bare values of the couplings. Let us consider the solution (|27|l in the neighbourhood 
of the UV cut-off A. Introducing a small positive parameter e = 1 — > 0, the term (j^)'^'^'^ can be expanded around 
e = 0: (1 — e)'^^2 = 1 — (d — 2)e + 0{e^). Inserting this into the exponent on the r.h.s. of Eq. 127|l and using 
ln(l - e) = -e + ©(e^), one finds 

exp(-ad/32(A)n2e + 0(e2)) ^ exp (^ad P^AW ~ e) + 0{e')^ ~ (1 _ e)". ^'(A)"^ _ (28) 



d / , X ad/3^(A)n^ + 0(e^) 



Then, the solution H27|) can be recast in the form 

7i„(fc)={t„(A)(^-j (^-j . (29) 

One may be tempted to conclude that the couplings Un{k) are relevant, marginal, and irrelevant depending on 
(P{A)n'^ — ^ < 0, = 0, and > 0, respectively. This distinction is, however, made according to the hare value of the 

coupling 0^{A). The latter goes to infinity in the UV regime (pushing the UV cut-off A to infinity), so that in fact 
all couplings Un{k) are irrelevant in the UV regime. As the scale k is decreased, a progressively larger number of 
couplings with increasing index n become relevant. For example, as Fig. ^(a) illustrates, the coupling, iti(fc) which is 
irrelevant in the UV regime, becomes relevant along the RG trajectories. In the cross-over region the UV scaling laws 
lose their validity due to the non-linear terms of higher order. Including the latter and following the RG trajectories 
through the critical scale k^ down to fc — *■ 0, one does not find any qualitative change in the behaviour of the RG 
trajectories. This is illustrated in Fign(b) using the example of a single RG trajectory (see Sect. IIII El for the details 
of the calculation). 

The UV scaling behaviour found for d > 3 is completely different from that obtained for d = 2 in Ref. ^2 [see 
Eq. (17) ihid]. For d = 2, a strong and a weak coupling phases with different scaling laws have been detected in the 
UV regime depending on the actual value of the scale-independent parameter (3. For d > 3 the same UV scaling laws 
have been found for arbitrary values of f}. Neither the UV scaling laws, nor the full RG-flow diagram indicate the 
existence of different phases in the lowest order of the gradient expansion. The disappearence of the topological phase 
transition for d > 3 of dimensions is of course related to the fact that long-range order is no longer prohibited by 
the Mermin- Wagner theorem in higher dimensions, and to the interplay of the long-range order with the topological 
excitations. 

E. Effective potential 

We now determine numerically the effective potential Ves{4>) for the 4-dimensional GSGM as the limit fc ^ of the 
blocked potential T4(0). For this purpose the flow equations (|21|) and (|23|l have to be solved numerically for k > kc 
above the scale kc of the spinodal instability, and the minimization steps according to Eq. H25|l have to be carried 
out for fc < fcc in the region of the spinodal instability. The Fourier-expansion of the periodic potential has to be 
truncated at the term with n — N, N fixed. Monitoring numerically the coupling constants w„(fc) (1 < n < N) of 
the periodic field theory, we compare the solution with the results obtained by solving the WH-RG equation lfT3|) in 
LPA for the appropriately chosen polynomial potential 

^ 1 

m=l ^ 

The polynomial potential was "matched" to the periodic one at the UV scale by equating the bare couplings g2m (A) 
to the corresponding coefficients in the Taylor expansion of the bare periodic potential at cf) — 0: 

N N r N 1^1 



(2m)! 



'-^^ (2m,) 
I— 1 ^ ' 



In order to avoid the spinodal instability for the polynomial case, the initial conditions for the polynomial potential 
were chosen specifically as = A) = 0.0001, g^^k = A) = 0.1 and gm^k = A) = 0, if m > 4. In this manner 
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for any given values of and N we get different sets of initial values m„(A) for the Fourier-amplitudes Un{k). For 
the comparison to the corresponding polynomial potential, the blocked periodic potential is again decomposed into 
a Taylor-series at each scale k. In this way, as a side-effect, one can investigate the influence of the violation of 
periodicity due to the Taylor-expansion on the RG flow of the periodic potential. 

Equations (|21l) and H23I) are solved numerically, starting from the UV cut-off A = 1 by using a fourth-order Runge- 
Kutta method with the step size 6k — lO^^fc with p — -i. There were no appreciable changes in the numerical results 
by increasing p further. The numerics have been tested by reproducing the results of '39] for the polynomial case. 
Below the critical scale kc where the spinodal instability occurs for the periodic case, the RG flow was determined by 
the help of the tree-level blocking relation (|25|l . As mentioned above, this equation has the same form for any number 
d of dimensions. Therefore, the numerical algorithm used and the function Pkiff) providing the minimum are the 
same as those for = 2 in our previous work 12]. Generally, for N > 0(10) Fourier modes a good convergence has 
been observed at all scales. For d — 4 the RG flow projected onto the plane (J, y) is qualitatively the same as that for 
d = 3 in Fig. ^ This flow diagram shows that no different phases occur in our approximation, but the UV- irrelevant, 
dimensionless parameter y becomes relevant in the IR regime. 

The upper plots in Figs. [21and|ni and those in Figs. 01and[3show the scaling of the dimensionful coupling constants 
52 (fc) and 54 (fc) above the scale of the spinodal instability for various truncations N of the Fourier expansion of the 
potential and for the bare values /3^(A) — 2/3q — 1287r^ and /3^(A) — ^/Sq — 327r^, respectively. [According to Eq. 
ui{k) is independent of the scale k for /3^(A) = $0 = 647r^.] We see the good convergence of the numerical results 
with increasing number TV of the included Fourier-modes. There is no significant difference in the extreme UV regime 
between the solutions for the periodic and the corresponding polynomial cases. Below this narrow UV "window" , the 
scaling laws for the coupling constants of the periodic potential deviate remarkably from those for the corresponding 
polynomial potential. Namely, the coupling constants for the periodic case start to decrease in magnitude after passing 
through a cross-over region. 

The lower plots in Figs. HJ— Elshow the scaling of the coupling constants 32 (fc) and (74 (fc) in the region of the spinodal 
instability for various values of the bare coupling /3^(A). The corresponding curves in the lower plots have not been 
artificially smoothened out in order to demonstrate the relatively small inaccuracy of the numerical calculation. In 
this deep IR regime the dimensionful coupling constants tend to zero for the limit fc — > in an oscillatory manner. 
Increasing the number N of the encountered Fourier-modes of the blocked potential does not alter this general 
behaviour. This behaviour is qualitatively different from that for the polynomial case when all dimensionful coupling 
constants (with even indices) tend to a non- vanishing constant for fc ^ |2^ [s^ • 

As illustrated in Fig. El the blocked periodic potential flattens out to a constant in the limit fc ^ 0. This behaviour 
is in accordance with the global features of the RG flow discussed in Sect. IIII Bl In this case the information on the 
IR physics of the system is encoded in the manner in which the blocked potential tends to a constant value. Because 
the spinodal instability always occurs, the minima of the blocked potential are filled due to the large-amplitude 
fluctuations of the field. 



IV. RG FLOW OF CG 



A. RS-RG for CG 



In this section we apply the RS-RG method in the dilute gas approximation (DGA) to the renormalization of the 
d-dimensional EGG. The same approach was applied to the two-dimensional generalized Goulomb gas many years ago 
in 0. Here, we generalize it for the the partition function (jSJ with in d > 2 dimensions while following as close as 
possible the derivation in 7] for d = 2. The details of the renormalization procedure are described in Add. IXI The 
RS-RG in DGA provides the following flow equations for the parameters of the d-dimensional EGG: 



d^ 
da 

dwa 
da 



-{d-m'+P'^[l2]<^">^ 

lad - ic2 5d,2 P'^ Wcr Cp- -\- ifi<J ^ \ 

a' ^(j 



exp 



-i^V'(a-a')G(a) 



-Wc 



Wa'W-ai exp 



i/3V'G(a) 



(32) 
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where 



2 \ ^ 2 

< a >= o" WaW-<y exp 



(33) 



and G{a) = (1 — 5d,2)[{d — 2)rirf]~^. Note that the first and second rows of the curly bracket on the r.h.s. of the first 
of Eqs. (|32|l correspond to the cases d = 2 and c? > 2, respectively. 

The well-known flow equations for d = 2 7] are recovered by these general formulae. For d = 2 and Z{2) symmetry, 
W-a = Wa, one gets G{a) — and < ct^ >= 2^^^ ct^w^ct, and the RG equations take the form: 



'^^^ O fl4Y^ 2 2 ( 2 

a—l ^ 



— Ait Wn 



<y' = l 



(34) 



Up to numerical prefactors, these equations are the same as Eq. ^ in Sect. Ill Bl which represent the well-known 
RG equations for the two-dimensional electric Coulomb gas 0- The extra factor 2 on the r.h.s. of the first equation 
appears due to the contributions of both charges cr = ±1 to the value of < cr^ >. The other, slightly different factors 
depend on the different inaccurate determinations of the overlapping regions of integrations for momenta, when the 
contributions of the annihilation and fusion of charges are calculated (e.g. see the note after Eq. ljA14p 'l. 



COMPARISON OF RG FLOW FOR RELATED MODELS 



RG flow for GSGM and EGG 



The fiow equations for the d-dimensional EGG can be compared to those for the GSGM. Such a comparison is 
meaningful only for the UV regime due to different renormalization schemes used in the two cases. Furthermore, the 
RS-RG method has the drawback that is not adequate to detect the spinodal instability. This is because the RS-RG 
applies a sharp cut-off in real space which corresponds to a smooth cut-off in momentum space. Therefore, there is 
no way to treat modes above and below the critical cut-off kc differently and to detect the spinodal instability in the 
RS-RG approach. 

In order to compare the RG fiow of the EGG with that of the GSGM, we rewrite Eqs. ^A2\ in terms of the 
dimensionless fugacities u„ — a'^u^ by making use of Eq. ©. We also introduce the momentum scale k cx a^^ . Then, 
the flow equations obtained by the RS-RG for the EGG take the form: 



k- 



" dk 

dUn 

Hk 



+ (d-2)/32- (1-3,5^,2) 



Ad 



a^UcrU-cr SXp 



(T=l 



^17 ^(7^(7 



exp 



'2^2 



/3V 



-(l"<5d,2) 



2(d- 2)f7d 
2{d~2)nd 



(35) 



~Ad 



■u„ 



Ua'U^^r exp 



-(1-^2) 



/32 a'' 



2{d~2)nd 



(36) 



These equations should be compared with Eqs. 1)21(1 and l(24|l for the GSGM. As to the tree-level scaling, i.e. the 
linearized fiow equations, the WH-RG LPA and RS-RG DGA give identical results. The non-linear term of the order 
0{u) on the r.h.s. of Eq. 1(24(1 is identical with the corresponding term on the r.h.s. of Eq. 1(36(1 for d = 2 and has 
a slightly different numerical prefactor for d > 3. The two methods give completely different results concerning the 
non-linear terms of the order higher than 0{u). This happens because the gradient expansion and the dilute gas 
approximation cannot be related directly. For comparison, the RG fiow may be calculated for both models in d = 3, 
and the RG trajectories may be projected onto the plane (J = 1//?^, y) where y = u\ for the GSGM [Fig. ^b)] and 
-0^ /{d>Ti)] for the EGG [Fig. [7]. In the latter case the projected RG trajectories are given by 



y = {2tt) 2{tf exp[ 
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A significant difference arises due to neglecting tfie wave-function renormalization in the WH-RG approacfi for the 
GSGM. Namely, Eq. (|21|l valid for the scale-dependence of the 'temperature' at any scale k, does not contain any 
loop-corrections as compared to the corresponding Eq. H35(l for the ECG. Therefore, Eq. H21(l does not have any 
fixed point at finite values of J and y (i.e. those of and ui) in the flow diagram for the GSGM [Fig. H^b)]. There 
exists a UV fixed point (J — > 0, y — > oo) and a zero-temperature fixed point (J — > oo,y —^ oo). The phase diagram 
for the ECG (Fig. [TJ however reveals the existence of a repulsive UV fixed point at finite values of J and y, when 
the terms of the order 0{y^) are neglected in Eqs. (|37b|) . For small values of the fugacities the non-linear terms are 
negligible and the RG flow looks rather similarly for both models. The fixed point in Fig. [7| maybe artificial because 
on the one hand, the RS-RG method in the dilute gas approximation is not well suited to determine the RG flow for 
large values of the fugacity. On the other hand, the WH-RG in the LPA, the validity of which is not restricted to 
small fugacity values, does not exhibit any flxed point. Nevertheless the clarification of the role of the wafefunction 
renormalization in the WH-RG method is an open question and it seems not to be possible in the framework of the 
gradient expansion. 

The comparison of the fiow diagram in Fig. Qwith that in Fig. |S1 where the terms of the order 0{y^) are also 
included, reveals the role of the last two nonlinear terms on the r.h.s. of Eq. I|37b|) . Due to these nonlinearities, 
the fixed point at finite parameter values turns from a repulsive UV fixed point to an attractive IR one. As to the 
global features of the flow diagram, this indicates that one certainly has to go beyond the DGA in order to clarify the 
essential role of the nonlinearities in the RS renormalization of ECG. 



B. RG flow for ECG and VLG 



The RG-flow diagram for the ECG presented above may be significantly modified due to non-linear terms of higher 
order when going beyond the DGA. Nevertheless, it can be compared in the UV regime to the RG-flow diagram for 
the VLG that was also obtained after certain linearizations |2Qj. Such a comparison could in principle give more 
insight to the question posed in the literature (see ^3 and the references therein) on whether the VLG and the 
three-dimensional XY-model (or the three-dimensional Coulomb gas) really belong to the same class of universality, 
or not. Our strategy is simply to show that there exists a transformation of the parameters {f3'^,Ui) i-^ {K,y) which 
transforms the set of Eqs. H35|l and H3()|l with the charges cr = ±1 for the three-dimensional ECG to the equations of 
the form very close to that of the set of Eqs. (|13|l for the VLG up to the non-linear terms of higher order. In fact, 
introducing the parameters 



the fiow equations give 

dK 



da 



-) 



^ K- 



47r^ 



1 

'4^ 



Ui exp 



Stt 



I— 

da 



= 6- 



K 
8^ 



(^)>.o(.^)^ 



(38) 



(39) 



These equations are very similar to Eq. H13|) . although there is an additional scale dependence in the equation for 
the flow of the parameter K. More questionable is the interpretation of y as fugacity of anything because it is a 
negative quantity. These discrepancies may reflect the dependence on the different approximations used when the 
RS-RG has been applied to ECG and VLG, respectively. Therefore, any conclusion on whether the VLG and the 
three-dimensional Coulomb gas belong to the same universality class, has to be seen as premature at this stage. 



VI. SUMMARY 



A renormalization-group analysis of the generalized sine-Gordon model (GSGM) and that of the equivalent Coulomb 
gas (ECG) have been pursued for d > 3 dimensions. The RG equations for the GSGM and the ECG are derived 
in two ways: (i) by means of Wegner's and Houghton's (WH-RG) method in the local-potential approximation and 
(ii) by the real-space RG (RS-RG) in the dilute-gas approximation. These different approaches complement one 
another: the WH-RG applied to the GSGM is better suited to determine the IR behaviour of the system, whereas 
the RS-RG applied to the ECG is able to better represent the non-linearities which influence the UV scaling. It 
has been shown that the leading-order terms of the flow equations obtained by these methods agree well. Taking 
as an example the 4-dimensional GSGM, we have demonstrated numerically that the blocked potential tends to a 
constant effective potential, independent of the field- variable, in the IR limit. This behaviour is in agreement with 
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the global requirements on the RG flow that constrain the effective potential to be periodic and convex at the same 
time. Comparing the evolution of the periodic blocked potential to that of a polynomial potential obtained by a 
Taylor-expansion of the bare periodic potential, we show that a violation of the periodicity would lead to a completely 
incorrect IR scaling. The flow diagram in Fig.^shows that only a single phase occurs in our approximation, and that 
the UV- irrelevant, dimcnsionless parameter y becomes relevant in the IR regime. 

Because the GSGM, the EGG and the XY-model belong to the same universality class, the determination of the 
RG flow of these models for the three-dimensional case via two different, complementary methods in the UV regime 
enables us to compare the UV scaling with that of a three-dimensional vortex gas (VLG) 20] . It has been established 
that one may transform the flow equations for the EGG to a form very close to that of the flow equations for the 
VLG, if only the leading order terms are kept. Nevertheless, the negative fugacity one needs to introduce and the 
different explicit scale-dependences occurring in the two sets of the flow equations do not enable one to make any 
definitive conclusion regarding the question whether the VLG and the three-dimensional Goulomb gas belong to the 
same universality class. 
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APPENDIX A: RS-RG FOR THE CG 



Here, we perform a real-space renormalization in the dilute gas approximation (DGA) of the CG given by the 
partition function ■ In order to follow as close as possible the procedure employed for d = 2 (see ^3 and references 
therein), we use continuum integrals over the infinitely large volume Dq, but cut out spheres Kj of radii a at each 
Coulomb charge aj centered at Xj (this is different from lattice regularization). Thus, the regions of integration for 
the positions of the charges are Dq, Di = Dq — Ki, D2 = Dq — Ki — K2, . . ., -D,y-i = -Do ^ Ki — ... — K^^^i for the 
charges cri, a2, ■ ■ (Ju, respectively. 

The RG flow equations can be found by changing the resolution, i.e. the smallest separation distance from a to 
a + (5a at which two charges can still be detected separately. This means that if two charges become closer than the 
separation distance a -\- 5a they will be seen as a single charge identical to the sum of the two original ones. The 
partition function will evolve due to the following reasons: 

1. the rescaling of the length, 

2. the fusion of charges in a single, non- vanishing one, 

3. the annihilation, i.e. fusion of charges in a single, neutral one. 

These sub-processes acting on the configurations with the number u of charges (described by Z^^^) give contributions 
5rsZ^'^\ 5fZ^'^~^\ and SaZ^'^^'^^ to the configurations with the numbers v-, (i^ — 1)-, and (j^ — 2) of charges, respectively, 
i.e. 



(Al) 



v=0 



u=2 



The RG idea is realized by recasting the modified partition sum Z + 5Z in the original form observed by Z , but with 
modified parameters 0^ and Wa- Below we determine all these contributions. 

1. Rescaling. Rescaling means the change of all integration variables Xj to x'^ — x.j{l + 5a/ a) and increasing the 
size of all regions Dj, Kj correspondingly. Then we rewrite the partition sum as an integral over the original 
variables and regions. First of all, the rescaling the volume elements df'x increases them by {6a/a) dd'^x; this 
can be merged into a modification of the fugacities, 



5rs,lWa 



5a 

Wad . 

a 



(A2) 



The change of the propagator should be treated differently for d = 2 and d > 3. 
The case d = 2. The propagator changes according to 



5rsG{x) 



5a 

a 



(independent of x). In the interaction energy this yields a contribution 



(A3) 



(A4) 



However, for neutral configurations we have X^J^i tTj — 0, i.e. J^i^j '^i — ~^3^ that one gets the energy 
contribution 



1 5a 



(A5) 



which leads to the following modification of the fugacities: 

5rs,2Wa = Wa exp --C2P CT 

\ 2 a 



1 fl2 2 
'7:C2P Waff . 

2 a 



(A6) 



The case d > 3. The change of the propagator, 



5^^G{x) = -{d-2)G{x)-^, 



(A7) 
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is proportional to the propagator itself. Therefore, it can be encoded in the change of the coupling 0^: 
Thus, we get for any numbers > 2 of dimensions 

drsP = -{d-2)P , drsWa- = Wad -C2()d,2P WaCr . 

a a 2 a 



(A8) 



(A9) 



2. Fusion. 

Let us denote the sphere of radius a + 5a centered at Xj by K'j and the annulus of inner radius a and thickness 
5a centered at Xj by 5j. Then the regions of intergration in the rescaled partition function are Dq, Di = 

Dq — K[ + 5i = D'l + 8i, D2 = -Do ^ K[ — K2 + 5i + 82, etc. for xi, X2, etc. respectively. The contributions of 
both the fusion and annihilation of the point charges ap and ag are included in the expression 



_ Y^V^ f d'^Xl r d'^Xq-l f d'^Xq r d'^Xq+l f d'^Xp^l 

r d'^Xp+i r d'^Xi, 

Jd'~^ Jd: 



d'^x 



Hence, we have to calculate integrals of the form 



f d'^x f d'^y , ^ 

JDo-K^ 0- J 6a: ^ 



(AlO) 



(All) 



where the point z satisfies either [.r — z| > 2a for /i, or a < |a; — z| < 2a for I2. Therefore, the integral for y in 
I2 extends only over the part 5x of the annulus 5x lying outside of K^: 



_ f d'^x r 

J\x-z\>2a Js.. 



d^x f d'^y 

" ■i''<\x-z\<2a " 



r d'^x f 

Ja<\x-z\<2a JSa, 



d'^y 



l\x-z\>2a 0- 

Introducing the center of mass coordinate t= h{x + y) and the separation coordinate s = x — y, one gets 



(A12) 



f dH^ 5a ^, 

\ —T^ld—j{t,s = a), 



dHl^ 5a ^, 
—^-Ud—f{t,s = a). 
a<|t|<|a a'' 2 a 



Then, the integral I can be estimated as 



f If 1 dH^ 5a ^, , f 5a ^, 



a). 



(A13) 



(A14) 



As to the regions of integration, we make another approximation: in accordance with the relation Kt{\a) C 
Kq U Kp C Kt{'^a)^ (the arguments indicate the radii of the spheres) we replace Kq U Kp by K^. Due to these 
approximations terms which are polynomials in a may be neglected. Therefore, our approximations will cease 
to work as o increases, so that our results are restricted to the UV scaling region. In order to take all terms of 
the order 0{^) into account, we perform a Taylor-expansion of the function f{t, s) in the variable s before the 
integration over the annulus. Thus, we find 



V p—l 



' — ' -. I. ^ ^ 



v=0 



d'^xi 



i-p=2 q=l''Do 

i 

Xp—l 



d'^x. 



q-l 



D' 



dH 



d'^Xq+l 



r d'^xp^i r 

Jd'^_, a'' Jd 



d'^ 



Xp+\ 



5a 
+ — 



E 



(7l,...,CTi,5^0 



Ja<\s\<a+Sa " \ ^ / 



. (A15) 
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In order to evaluate this, we expand the exponent first, 



T<Ti ,...,(T„ 

Xl,...,Xiy 



j[""'lcxp|-i/3Vpa,G(a)|exp|-;32 ^ ^^^^ 



1 - \so.dl + ^SaS/3a^9^ 



1 + + ^saspdid}, 



(A16) 



Expressions (|A15|) and (|A16|) represent the starting point for the determination of the contribution of the fusion 
and annihilation of charges to the partition sum. In case of fusion, when the charges ap and aq fuse into a single, 
non- vanishing charge fft ~ <^p + '^q =/= 0, the only contribution of the order 0{5a/a) comes from the zeroth power 
of s in HA16() . Renumerating the charges in the manner that the new charge at i = x^^i becomes the {v — l)-th 
one, we can write 



D' 



a< I s| <a+(5a 



Sa 



D' 



(A17) 



The other terms with Jg s"^d'^s ^ 0(a™)^ where m — 1,2, .. . are neglected. Inserting this in the change of 
the partition function due to fusion. 



(5a 



-l]g-i/3Vp£T,G(a) 



d'^Xl 



A calculation of this expression finally yields 



^f^o- = ^^"i"" E '^<y'^<y-<y' exp 



(A19) 



Thus, fusion generates a renormalization of the fugacities. In the manipulations above, we considered all charge 
configurations as identical if they can be obtained from a given one via permutations only. 

3. Annihilation. 



The charges ap and Oq ~ —Up fuse into a neutral particle when the inequality a < \xp — Xq\ < a + 5a holds for 
their distance. This is called annihilation, since both charges vanish from the gas. Effectively, the term Z^'^^ then 
yields a contribution to Z^'^~'^\ In Eqs. ljA15|l and (|A16p one has to keep the terms quadratic in the separation 
coordinate s (those of the order a?'^'^^'^^ in the exponent of jl'^l) in order to find all terms of the order 0{5a/a). 
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Instead of the integral ljA17|) , we now get 



D' 



^'^ J a<\s\<a+5a ^'^ 



^ J a<\s\<a^5a ^ 



Xj . 



j[u-2] 



exp 



1 



-/3V^G(a) 



D:,_^ ^'^ Ja<\s\<a+Sa '^'^ 
1 



(A20) 



Here the hnear term and the off-diagonal part of the matrix SaSfj give vanishing contributions to the integral 
over the annulus, so that one gets 



dH 



d'^s 



D' „ J a<\s\<a+5a " 



jl''-^] exp 



1 % 



dH r d<^s 

^1-2 " Ja<\s\<a+Sa O, 



(A21) 



We now change the order of integration so that the annihilating charges are the "last ones in the queue" and 
integrate by parts in the second term, with the result 



dH 



" J a<\s\<a+Sa 



exp 



1 



5a 



iy^2 



1 



/3'«'"'^P E J^, (A22) 



where the last integral vanishes since ddG{x) ~ S{x) and the point t — Xi is excluded from the region of 
integration. The first term on the r.h.s. can be included in the renormalization of the free energy which, 
however, we do not consider here. Therefore, the only nonvanishing terms originate from the surface integrals 
on the r.h.s. of Eq. (|A22() . They should be considered separately. The integral over the surface of the domain 
Do vanishes as the integrand goes sufficiently fast to zero for \t\ — > oo. The integrand of the integral over the 
surface dKi of the sphere centered at xi contains the derivative 



(t - Xi)c 



-C2, 



for d=2. 



\t-x^\'^ I i(d-2)G^^,t, for d>3. 
Using the outer normal tIq, of the surface dKi and writing t — xi + an, we obtain: 



dKi 



d(JaGxj,t d^Gxi 



df naG{xi + an — xj) d\^G{xi + an — Xi) —: I . 



(A23) 



(A24) 



dK, 



Here we have to distinguish the following particular cases: 
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(a) For / ^ ij: 



dfua 



dKi 



G{xi - Xj)dlG{xi - Xi) 



+ani3d'i^G{xi - Xj)d^^G{xi - Xi) + anpG{xi - Xj)djjdlG{xi - Xi) 



(A25) 



Here, the first term vanishes. The last term vanishes too as it is proportional to 5{xi — Xi), but we have 
explicitly xi ^ Xi. The second term can be evaluated as 



^i—d\G{xi - Xj)d'^G{xi - X,) = rid 



Q*^ {xi,j)a{xi^i)a f c|, for d = 2, 

d koP^Mp U(rf-2)2G(xi,,)G(x/,,), for d>3, 



(A26) 



with xij = xi — Xj. 
(b) For l=j^ i: 



d<TaGxj,t dlGx^^t 



df Ha G{an)d]^G{xi + an — Xi) 

dKi JdKi 

I ril — N I , 



dKi 



dfuo. G{an)[d'^Gixi,,) + anp d'^d'^^Gixi,,)] = 0, 



(A27) 



where the first term vanishes due to isotropy and the second one vanishes, too, for the same reason as the 
last term on the r.h.s. of Eq. ljA25|) . 

(c) For I = i ^ j: 



I ^ f daaGa:^,tdlGx^.f 
IdKi 



dfria G{xi +an~ Xj)d^G{y) 

dKi V=a-n 



dKi 



dfn^ G{xi,j) + anpd'f,G{xi,j) dlG{y) 



dKi 



for d = 2. 



{d~2)G{an), for d>3. 



C2 , for d = 2 , 

2(d-2)G(a), for d>3. 



(A28) 



(d) For I = i = j: 



d(JaGxi,tdt,Gxi,t 



dK, 



dfnaG{an) dlG{y) 

dK, !/=an 



0, 



for d = 2, 



^Qda'^-^id-2)G^{a), for d>3. 



(A29) 
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By employing all the above results, we are now finally in a position to rewrite the r.h.s. of Eq. ljA22|l 



1 s 



J exp 

v-2 v-2 



fid 



Sa 1 
a a'^ 

v-2 



TV - (v-2) 



j,i=l l=^ij 



J" 



C2 

i(d-2)G(a) 



4 

!/-2 



(A30) 



where, again, the first and second rows of the columns stand for d — 2 and cJ > 3, respectively. Multiplying this 
expression by WapW-a,, and summing up for cTp, one obtains: 







a 






p 1/-2 


v-2 V 


-2 




3 1^3 



WcrW-o- exp 



v-2 



l^^,3 



1 (Xij)a(a;;,i) 

a I O' 



C2 

2)G(a) 



v-2 

Y,^h^da-\d~2)G\c 



(A31) 



with 



< >~ a'^Wa-W-cr exp 



i^V2G(a) 



(A32) 



With the help of this relation, we can now write down the infinitesimal additional contribution to the partition 
function under the blocking transformation. 



Z + daZ = ( 1 + rid— ^-7- WcrW-a- exp 



d'^xi 



d 



Dr a 

d — 2 Sa 
a 



d'^x^^2 



\i3^a^G{a) 



E 



1 v\ 



v=2 ^ ' 



D',_ 



n E '^-^ ( ^ " E ^'^p 



■j = l (Tj 



\p'a^G{a) 



~A Sa ( -C2 



^ ^ d a \ i(d-2)G(a) J 2 



!/-2 1^-2 



3 1^3 



(A33) 



Here, the first term leads to the renormalization of the free energy. The contribution of the third term on the 
r.h.s. of Eq. ljA31|l vanishes that can be seen after partial integration since xi^i and xi^j are not zero. The effects 
of all the other terms merge in the modifications of the parameters of the original partition function: 



SaP'' 




d a I 




Sa n^. 




= Wa-j- 




a d 



-C2 



Wa-'W^cr' exp 



i/3V''G(a) 



^P^wJ-nla^GHa) <a'^ > . (A34a) 
4d a 



As we sum the various contributions to the changes of the parameters, we finally arrive at the fiow equation H32|l . 
which completes its derivation. 
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J J 

(a) (b) 



FIG. 1: RG flow for the three-dimensional GSGM, obtained by solving the 
WH-RG equations (a) Eqs. and (b) Eqs. ^ outside the 

region of spinodal instability (solid line) and Eqs. 121L I25I I inside the re- 
gion of the spinodal instability (dotted line) by a fourth-order Runge-Kutta 
method. The RG trajectories are projected onto the subspace of the dimen- 
sionless coupling constants y = Mi(fe) and J = l/f3^{k). The lines represent 
the RG trajectories for various initial conditions, and the arrows indicate the 
direction of the RG flow. 
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FIG. 2: RG-evolution of the dimensionful coupling constant 32 (fc) versus the 
running momentum cut-off k (from k = 1 ... 0) in d = 4 [as a starting point 
for the RG flow, we fix /3^(A = 1) to the value 1287r^]. Various truncations N 
of the Fourier-expansion of the periodic potential have been employed. The 
upper and lower plots show the scale-dependence above and below the scale 
fee of the spinodal instability which occurs at kc — 0.0022 for A'^ — 10. The 
solution for the corresponding polynomial case is also indicated (solid line) 
MM. 
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FIG. 3: The same as Fig. |21for the coupling constant 34 (fc). 
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FIG. 4: The dimensionful coupling constant 32 (fe) is plotted versus the run- 
ning momentum cut-off k (from fc = 1 to 0) for d = 4 dimensions, with 
/3^(A = 1) = 327r'^. Various truncations A*' of the Fourier-expansion of the 
periodic potential are employed. The upper and lower plots show, respectively 
the scale-dependence above and below the scale fcc of the spinodal instability 
which occurs at fee = 0.01413 for N = 10. The solution for the corresponding 
polynomial case is also indicated (solid line) [2^l39| . 
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FIG. 5: The same as Fig. |l]for the coupling constant g4,{k). 
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FIG. 6: The flattening of the dimensionful periodic potential Vk{4>) is plot- 
ted inside the region of spinodal instability, k < kc- The various curves with 
decreasing curvature at tf) — correspond to different values of the scale de- 
creasing from fe = fee in steps of Sk = fcc/10. 




FIG. 7: RG flow for the EGG in the plane (J, y), calculated without the terms 
of O(y^) on the r.h.s. of Eq. ll37bl . 




